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We consider how Higgs Inflation can be observationally distinguished from variants based on gauge 
singlet scalar extensions of the Standard Model, in particular where the inflaton is a non-minimally 
coupled gauge singlet scalar (S-inflation). We show that radiative corrections generally cause the 
spectral index n to decrease relative to the classical value as the Higgs mass is increased if the Higgs 
boson is the inflaton, whereas n increases with increasing Higgs mass if the inflaton is a gauge singlet 
scalar. The accuracy to which n can be calculated in these models depends on how precisely the 
reheating temperature can be determined. The number of Einstein frame e-foldings A'^ is similar in 
both models, with iV ~ 58 — 61 for singlet inflation compared with N 57 — 60 for Higgs inflation. 
This allows the spectral index to be calculated to an accuracy An — ±0.001. Provided the Higgs 
mass is above « 135 GeV, a combination of a Higgs mass measurement and a precise determination 
of n will enable Higgs Inflation and S'-inflation to be distinguished. 

PACS numbers: 



I. INTRODUCTION 



At present our knowledge of particle physics is entirely summarised by the Standard Model (SM). There is no direct 
evidence of particle physics at mass scales greater than the weak scale, although the SM leaves many unanswered 
7*^' questions, such as the origin of the baryon asymmetry, the nature of dark matter, the solution to the strong CP 
Qh' problem, the detailed structure of gauge and Yukawa coupling constants, the origin of neutrino masses and a particle- 
^ [ physics interpretation of inflation in the early universe. Nevertheless, since the only particle physics mass scale we 
know to exist is the weak scale, it makes sense to explore the possibility that the answer to these questions can be 
obtained from models of weak scale particle physics. This has the added appeal that such models may be testable at 
colliders. 

' Two of these questions might be solved by scalar fields with masses and interactions characterised by the weak 
scale. Firstly, Higgs Inflation [1] proposes that the SM Higgs boson could be responsible for inflation. This is possible 
if there is a large dimensionless coupling of the Higgs to the Ricci scalar through the term £^hH^HR, with ^ 10^. 
[ Secondly, a real gauge singlet scalar interacting with the SM via the 'Higgs portal' provides a minimal extension of 
the SM with a thermal relic WIMP dark matter candidate 0-01 • In this case there is a second possibility for inflation, 
since the singlet scalar could also have a non-minimal coupling to gravity and so provide an alternative to Higgs 
Inflation, which we call 5- inflation [t^. Higgs Inflation in a gauge singlet extension of the Standard Model has been 
discussed in Q. 

In all of these models the scalar couplings are O(O.l) and the potential is made flat by choosing a large non- minimal 
coupling S^fh of the inflaton^ (j) to gravity, such that 

X 



^ 3(0.027)4 



where is the quartic self-coupling of the inflaton and N is the number of Einstein frame e-foldings between the 
WMAP pivot scale leaving the horizon and the end of inflation. This large non-minimal coupling to gravity causes 
the theory to be effectively scale- invariant at large values of ^ « O(0.1)Mp. Thus slow roll inflation can take place. 
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^ In the following, is used to refer to the inflaton in general, which can be either s or h, depending which version of the model we are 
considering. 
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It is remarkable that inflation models based on weak scale particle theories have the predictive power to connect 
observables from cosmology (the spectral index) and particle physics (the Higgs mass). It is therefore important to 
clarify the predictions for the spectral index as a function of the Higgs mass ruh, as the models discussed in this paper 
will be within the reach of Planck and the LHC in the next few years. 

At the classical level, Higgs Inflation and its singlet variants make the same predictions for the inflation observables. 
At the quantum level, the spectral index of the models as a function of Higgs mass will differ. The predictions for 
Higgs inflation have been discussed in 0-13 • Our goal in this paper is to calculate the predictions for the spectral 
index and its dependence on the Higgs mass in the different models using the same method throughout. In this way 
we can be sure that any differences in the predictions are due to the fundamental differences between the models 
and not due to the initial conditions, method of obtaining the effective potential or specifics of the code used. We 
will show that a combination of Higgs mass measurement and a precise determination of n can distinguish between 
S'-inflation and models where the inflaton is the Higgs boson. 

It has been suggested that Higg s Inflation and its singlet variants may not be consistent particle physics theories 
at the quantum level (see |15l - [20l | and references within for discussion). It is known that perturbative unitarity is 
violated in tree-level graviton exchange-mediated scalar particle scattering via the non-minimal coupling to gravity^ 
at _E « A = ^ [l^. This scale A is below the scale of inflation. Therefore, something must happen at the scale A 
— either strong coupling ensures that true non-perturbative unitarity is not violated even though it is at tree-level, 
or some UV completion of the theory must restore unitarity. If new physics is necessary, then the model is no longer 
natural as we would need to add terms to the Lagrangian which would be large at the scale of inflation. However, 
if the theory is non-perturbative with no new physics, then we can make predictions from inflation, as this is not 
directly related to particle scattering processes. The possibility of strong coupling as a solution to unitarity has been 
discussed in [ill HJ^. 

Therefore, if strong coupling unitarizes the theory, Higgs Inflation and its variants would be consistent theories, 
requiring no new particles or interactions*. This can be tested by the resulting predictions of the model. Our 
philosophy for the remainder of this paper is to consider the possibility that Higgs inflation and its variants are 
consistent theories and therefore should be subject to rigourous testing against experimental data. 

Our paper is organised as follows. In Section 2 we introduce the three models of non-minimally coupled inflation 
that we consider: (i) original Higgs Inflation with no additional particles beyond the Standard Model (ii) inflation 
along the Higgs direction in an extension with a gauge singlet scalar Q and (iii) S'-inflation Q. In Section 3 we 
compute the quantum effective potential for the models. In Section 4 we discuss reheating. This is crucial, as the 
ability to accurately predict the reheating temperature is a key feature of this class of models, allowing the number 
of e-foldings of inflation for a given length scale to be precisely determined and so the spectral index to be precisely 
predicted. In Section 5 we compare the predictions of Higgs Inflation and S'-inflation at the quantum level. In Section 
6 we present our conclusions. 

II. NON-MINIMALLY COUPLED MODELS OF INFLATION 

Higgs inflation action in the Jordan frame 

We consider three variants of the Higgs inflation model, which can all be described by the following action for 
different choices of the couplings and inflaton field. The Jordan frame action, including quantum corrections, is 



^ It has also been suggested that the non-polynomial potential in the Einstein frame could also lead to unitarity violation at E A. We 
believe that this will not be the case if the non-polynomial potential of the inflaton interpolates between renormalizable polynomial 
potentials, since we would expect scattering processes at i? ^ A be dominated by the renormalizable potential at large </> ~ £ up to 
small corrections proportional to powers of A/E, vanishing as A — >■ 0. 

^ An argument against this point of view is presented in |23 . |23|I . 

* However, if unitarity is really violated in scalar scattering processes, it may still be possible to add interactions and particles to unitarize 
the theory HIl^. 
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where C-^^ is the Standard Model Lagrangian density minus the purely Higgs doublet terms and V{H^H, s^) is the 
renormalisation group (RG) improved effective potential, given by 

V^X, [g% (H^H) - + ^Gj^s^H^H +^ + (3) 

Xhs is the Higgs portal coupling. Gnit) — exp ^— /J and the anomalous dimension is given by Eq. (jA-7|) . 

while = 1 (as 7s =0). For 'pure' Higgs inflation, we set s = Xs = Xhs = = 0. We also consider the case 
where the inflaton is the Higgs boson in the presence of a singlet scalar setting ^s(?7it) = 0, and the case where 
the singlet is the inflaton [1], setting £,h{'mt) = 0. Setting the non-inflationary non-minimal coupling to zero is a 
reasonable simpliflcation. If and £_h are of similar magnitudes, then the potential in the Einstein frame (oc S,^^) 
would be similar in all directions and inflation would be expected to occur in a general direction in the s — h plane. 
Our aim is to compare the two limiting cases where h and s are the inflaton respectively, so flxing the non-inflationary 
^ff, to be small (and so the potential in that direction to be large), ensures that inflation will occur either along the 
Higgs or the s direction. 



Conformal transformation to the Einstein frame 



We will transform the entire RG improved action to the Einstein frame, where the flelds are minimally coupled to 
gravity. This will allow us to use the standard slow roll formalism to calculate the spectral index n and tensor-to- 
scalar ratio r. Quantities in the Einstein frame will be denoted by a tilde (e.g. 5^1/). For now we consider only the 

physical Higgs field h, where H =^ 1^ ^ real. The inclusion of the non-physical components of the 

Higgs field is important when computing quantum corrections to the effective potential and also in the analysis of 
unitarity-conservation, but they are not important for the dynamics of the infiaton field and the calculation of the 
spectral index. 

The Jordan frame and the Einstein frame are related by a conformal transformation which transforms the metric 
in a field dependent way. Considering the action Eq. for general h and s, the conformal transformation to the 
Einstein frame is defined by 



f]2 



(4) 



with 



p 



M2 

p 



(5) 



Then 
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where R is the Ricci scalar with respect to g^i^ and 
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dt 



(7) 



We can then rescale the fields using^ 



dXh 
dh 



ds 



(8) 



^ These are only total derivatives in the limits s — > and h respectively. 
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to give 



where 



and 



d^x^g(c^- + -r^'d^Xhd.Xh + ^r"" d^XsduXs + Aixs^XhYg^" d^^xud^Xs - U{xs,Xh)), (9) 



^(x.x.) = ^t^^^/^« (10) 



Uixs, Xh) - ^ {^^G\,h' + + ^GWh') (11) 

During inflation, only one field (the inflaton) is non-zero and so the non-canonically normalised term A is zero. After 

inflation, the Jordan and Einstein frames will be indistinguishable (as both <^ Mp and |- <C Mp). Therefore 
the curvature perturbation spectrum calculated in the Einstein frame can be compared to measurements made in the 
physical Jordan frame. 



Classical predictions from slow-roll inflation 

Using the tree- level potential and the approximation ^ 1, the tree- level slow-roll parameters are 

p 

4 M'' 4 1\/P 

g~ _ P ■ fi _ P (12) 

where (/)^ ~ AMpN /ZS^^ is the field at N e- foldings. A calculation of the classical spectral index and tensor-to-scalar 
ratio then gives 

rici^l-^ ^ -HO (^) = 0.965 ; (13) 



^ + o(^-^] =3.6x10-^, (14) 



where we have used N = 58. This is close to the WMAP pivot scale in these models. The tree-level potential and 
classical predictions are identical for all non-minimally coupled scalar field models with a (jpR coupling to gravity and 
(f)'^ potential at large (f). 



III. QUANTUM CORRECTIONS TO THE SCALAR POTENTIAL 

In this section we compute the RG improved effective potential for Higgs Inflation and S-inflation. We calculate 
the effective potential in the Jordan frame using the Fcyiiman rules for the SM fields defined in the Jordan frame. In 
this approach, the effect of the non-minimal coupling to gravity is taken into account by using a modified propagator 
for the non-minimally coupled inflaton field. This is derived by quantising the inflaton field in the Einstein frame but 
without rescaling the field to the canonically normalised form. (We review this in Appendix A.) Using these Feynman 
rules, the RG equations for the couplings are derived and the effective potential is calculated in the Jordan frame. 
This is then transformed to the Einstein frame, where it is used to study slow-roll inflation. While it is also possible 
to use the 1-loop Coleman- Weinberg potential, this would not be the standard flat-space form but must also include 
the contribution of the non-minimal coupling^ of </> to R. 



In our previous paper on S'-inflation, the contribution of the non-minimal coupling to the Coleman- Weinberg potential was not included, 
but this has a small eflfect in the case of a real gauge singlet scalar inflaton. 
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Initial conditions 



We choose £d ("if) such that the model is correctly normalised to the WMAP 7-year mean value for the curvature 
perturbation [23, , which requires 

J = (0.00275Mp)'^. (15) 

The initial values of the coupling constants are defined at the renormalisation scale fi = mt, with mt — 171.0 GeV 
and V = 246.22 GeV. The gauge couplings are given by 

- 0.03344, ^ 0.01027 and ^^0.1071. (16) 



The couplings g and g' are obtained by an RG fiow from their values at fi — Mz, which are given in (29| . while 53 is 
calculated numerically. (See [13] and references within for details.) The pole mass matching scheme is used to set the 
initial conditions Xh{mt) and yt{mt). This relates the physical pole masses to the couplings in the MS renormalisation 
scheme through the following expressions: 

2 

V2 

yt{mt) = — mt(l + AO (17) 

V 

where A/j and Af account for radiative corrections and are given in the appendix of (26j . 

For both methods we fix TV = 58 for the WMAP pivot scale (we will discuss the value of N in the next section), 
and use this to determine 0^ via 

U_ fdxi 



N=l Hdt^l ^^^^)#, (18) 



RG equations and slow-roll parameters 

The RG equations for yt, g, g' , Xh, As, Xhs, ^s, £,h and the commutator suppression factors Ch and Cs are presented 
in Appendix A. The suppression factors are inserted for each infiaton propagator in a loop. For (f> <C — 1 while 

for (f) 3> "c^, ca — )■ 0. Using the RG improved potential we calculate the slow roll parameters e and 77 analytically as a 
function of the running couplings, as discussed in Appendix A. The spectral index is given by n = 1 — 6e+2fi. We find 
that the spectral index is dominated by fj, which has only small radiative corrections. However, the tree level value 
of 77 depends on 0^, which is determined by e. The radiative corrections to e can be large and therefore radiative 
corrections can have a large effect on the spectral index. The parameter e is given by 



f d(l)_V fdUl_ 
Ml fdl_y f A 

2 [dxj Un^ <P 



(19) 
(20) 



In the second expression is the tree-level contribution and gives the approximate contribution from radiative 
corrections. During infiation (i.e. = 0), the expressions for are defined as 

and 



IQn'U ^ — ^ -Qyt + ^j- 6A„ - 2A,,g (21) 



IGn'L, ^ Ms _ f 12A. + 6y^ - ^ (g^ + g'^) ] (22) 



6 V 2 
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where the terms oc £,h/£,s and £,s/£,h are hkely to be subdominant. We will use the approximation Eq. (|20p to explain 
the main features of the results. However, the numerical results are obtained using the full expressions for n, e and ry, 
given in Appendix A. 

IV. REHEATING 

In most inflation models, the uncertainty in the reheating temperature introduces a large uncertainty in the value 
of N corresponding to the scales observed in CMB experiments and so in the theoretical value of the spectral index. A 
remarkable feature of the models we are considering is that the reheating temperature Tr, defined to be the effective 
temperature at which the inflaton energy density is equal to the energy density in relativistic particles^, can be 
determined to a high degree of precision, allowing a precise prediction of the spectral index. 

The mechanism for reheating in these models will give us a narrow range of values for Tr. This will enable us to 
calculate N to high precision - a necessary ingredient in calculating the spectral index n (Eq. ([T^ ). Unlike most 
inflation models, the couplings of the particles to the Standard Model in Higgs Inflation and S'-inflation are either 
known or well-constrained, making this calculation in principle possible to arbitrarily high precision. 

However, in practice, the calculation is difficult, as non-perturbative effects are important. Reheating in Higgs 
Inflation has been studied by two groups P, |3l|. There has been no similar study of reheating in S*- inflation (although 
some discussion was presented in [32|). We will show that the process is very similar to reheating in Higgs Inflation, 
with stochastic resonance of Higgs bosons to gauge bosons being replaced by stochastic resonance of gauge singlet 
scalars to Higgs bosons via the Higgs portal coupling. 

A concern in the case of ^-inflation is that the dark matter scalars s are stable and so reheating via inflaton decay 
is not possible. We must therefore demonstrate that reheating can occur and that there is no dangerous density of 
stable s scalars remaining from inflation. We will show that this is the case; reheating occurs via stochastic resonance 
rather than decay, and any remaining energy density in s oscillations will be thermalised by the dominant radiation 
background from reheating. 

We first summarise our results. Following the approximate analytical method of 9], we find that Tji is high for 
5'-inflation: 

3 X 10" GeV <Tr<8x lO" GeV (23) 
and similar to Tr for pure Higgs inflation [9| 

3 X 10" GeV <Tr< 1.5 x lO" GeV. (24) 

The small range of Tr for both models enables us to reasonably estimate N. For 5- inflation this is 58 < TV ;$ 61. In 
this we have conservatively included an additional error AN = ±1 to account for possible additional uncertainties in 
the theoretical estimate of Tu, giving a classical spectral index 0.965 ^ nd ^ 0.967. A measurement of Xhs through 
dark matter detection experiments would further increase the predictiveness of this model. For pure Higgs inflation, 
57 ^ TV ^ 60 (also including an additional theoretical error A7V = ±1), which gives 0.964 < rtd ^ 0.966. For 
the case of inflation in the Higgs direction including a singlet scalar, the result for Tr will be similar to pure Higgs 
inflation, although a little higher (as the inflaton has an extra channel to annihilate or decay into). We next review 
the mechanism for reheating in Higgs Inflation before discussing the process in S'-inflation reheating. 

Reheating in Higgs inflation 

An analytical calculation of was carried out in Q, while in [sTj a numerical calculation was performed. Back- 
reaction was discussed in 'sT] , but was not considered in since the annihilation of the gauge bosons produced via 
the resonance was found to be dominant in j9] . We now further describe reheating in Higgs inflation according to @ , 
as we will follow this method to calculate Tr for S'-inflation. 

Reheating in Higgs inflation occurs through a stochastic resonance process fij. After inflation, the Higgs-inflaton 
Xh oscillates in a quadratic (Xh) potential. The gauge boson masses are proportional to the modulus of the oscillating 
inflaton field. When the oscillation modulus is small during the oscillation cycle, non-relativistic gauge bosons are 



^ The temperature at which thermal equihbrium occurs may be lower than this |3C| . 
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produced non-adiabatically; when the modulus is large, these gauge bosons rapidly decay to relativistic Standard 
Model particles [9|]. This prevents the build up of gauge bosons and backreaction. However, the expansion of the 
Universe causes the maximum amplitude of the Higgs oscillation to decrease, which decreases the gauge boson mass. 
As the decay rate of the gauge bosons is proportional to their mass, eventually the gauge bosons can no longer 
decay appreciably and their occupation number builds up, enabling the stochastic resonance. At this point, the 
energy of the inflaton is quickly transferred to the gauge bosons, which in turn quickly annihilate to relativistic 
fermions. Reheating via the resonant production of Higgs excitations is also possible, but the stochastic resonance 
to gauge bosons is expected to dominate [9]. In Higgs Inflation, the resulting temperature of radiation domination is 
3 X 10^^ GeV < Tfj < 1.5 x 10^'' GeV (where the lower limit is from reheating via the production of excitations of 
the Higgs) 

The WMAP pivot scale is ko — ^ — 0.002 Mpc^^. The number of e- foldings at which this exits the horizon is 
given by 

where g{Tii) ~ 106.75 for the SM, g{To) ~ 2 and Tq is the present photon temperature. This assumes that the 
universe is matter dominated by inflaton oscillations from the end of inflation until the moment of reheating. For the 
range of from Higgs inflation, and including an additional theoretical uncertainty iV = ±1, we obtain 57 < TV < 60 
for the number if e-foldings at which the pivot scale exits the horizon. 

From the expression for the spectral index Eq. p3p . an uncertainty AiV corresponds to an uncertainty An = 
zLAN / N'^ . With AN = 4 and iV = 60 this gives An — ±0.001. Therefore, even with a conservative estimate of the 
uncertainty in N, the spectral index in Higgs inflation can be calculated to an accuracy ±0.001. This can be improved 
simply by improving the accuracy of the reheating temperature calculation; in principle, there are no obstacles to an 
arbitrarily precise calculation of the spectral index in these models. 

Reheating in S-inflation 

Reheating in S'-inflation occurs by the production of Higgs bosons via the coupling of s to H, which can then decay 
to the particles of the Standard Model. Since s is stable, any leftover inflaton density can only be transferred to 
thermal radiation through scattering with Higgs bosons and s bosons in the thermal background. We therefore need 
to check that the inflaton is a subdominant component of the Universe when its oscillation is small enough that it 
effectively oscillates in a quartic potential (oc xt ~ ■s*)- I^i this case, any residual oscillating inflaton density will 
eventually be thermalised by scattering from the dominant thermal background. Taking a conservative approach, 
we therefore require that reheating should occur during the early quadratic (xi) potential stage, where Xs is the 
canonically normalised Higgs field. (This is conservative as it still may be possible to have reheating once Xs is 
oscillating in the quartic potential.) Reheating occurs via stochastic resonance to Higgs bosons through ^^s^|-ffp. 
The process is very similar to the case of reheating in Higgs Inflation, with the gauge boson final state particles 
replaced by Higgs bosons and the gauge couplings replaced by Xhs- We therefore apply same method as used in the 
case of Higgs Infiation [9] to estimate the reheating temperature. Our calculation is outlined in Appendix B. 

The key points are as follows. Once the Xs oscillation amplitude becomes small enough to allow stochastic resonance 
to occur, the energy in the Xs oscillations rapidly transfers to the Higgs bosons, which in turn decay to SM fermions 
(primarily top quarks, since the Xs oscillation amplitude gives the Higgs bosons a large mass). This causes the Xs 
amplitude to rapidly decrease, such that 17 — > 1 and the s field becomes canonically normalised in an potential. 
Provided that the radiation energy density from Higgs boson decay is dominant at this time, the s oscillations and 
radiation will subsequently decrease as and the radiation energy density will therefore remain dominant at least 
until the time at which the s mass term comes to dominate the s oscillations. Provided that the thermal radiation 
background can thermalise the subdominant s oscillations before this time, there will be no stable s density remaining 
from inflation. In this case s dark matter will be a purely thermal relic density. 

In fact, there are two possible mechanisms for reheating in S'-inflation. The first is through stochastic resonance to 
Higgs bosons, which gives 

/A \ 

Tfl « 3 X 10" i^-^j GeV . (26) 
This requires As > 0.25Aft,s in order for reheating to complete before the potential changes from quadratic to quartic. 
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10-" 10-2 -L 

X s 

FIG. 1: Showing the regions of the A3 and Ajj^ parameter space allowed by the constraints of reheating. The shaded region is excluded 
under the conservative assumptions described in the text. In the majority of the region marked 'both', reheating via stochastic resonance 

is expected to occur first and therefore to dominate. 

Imposing a bound^ y^- ^ 10^ gives 

hs 

3 X 10^3 GeV <Tr<8x 10^^ GeV. (27) 

The lower bound is similar to that for Higgs inflation 9] . The upper bound for 5-inflation is higher because of the 
freedom in and Xhs- 

The second mechanism is reheating via the production of excitations of the inflaton itself, which subsequently 
annihilate to Higgs bosons. This gives 

Tfl « 9 X lO^^X]/^ GeV (28) 

and requires As > 0.019 to complete before the potential becomes quartic. 

A region of the parameter space is excluded due to the conditions on the couplings, as shown in Fig. ([T]). However, 
we are rather conservative in our estimates since some reheating is still likely to occur once the oscillations become 
quartic. With Tfj given by Eq. (I?f)) we find 59 < iV < 60. Including a conservative estimated theoretical error of ±1 
on N we arrive at 

58 < iV < 61. (29) 

Thus we expect the uncertainty in N and the resulting error in the theoretical estimate of the spectral index to be 
similar to the case of Higgs Inflation, An — ±0.001, which could be improved with more careful calculations. 

Thermalisation of the remaining inflaton oscillations 

Because s cannot decay, it is necessary to check whether there is any energy density remaining in the inflaton 
oscillations after reheating. From the end of inflation (when resonant reheating occurs) until the time when the 
mass term in the s potential comes to dominate the term, the radiation background will dominate the inflaton 
oscillations. The SM radiation background will thermalise well before the term comes to dominate. To check that 
the remaining energy in the s oscillations is thermalised, we can simply check that the energy in the oscillations, 
equivalent to a density of zero momentum s particles, is thermalised before the oscillations dominate the thermal 
background. 

The condition for sufficiently rapid thermalisation of the inflaton oscillations in the potential is F > 7J, where 
F is the scattering rate of zero-momentum s particles from particles in the thermal background (which also includes 



There are constraints due to vacuum stability and perturbativity of the couplings (see Q) and from the requirement that A^^ must be 
large enough to achieve the correct relic density of dark matter. 
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thermal s particles) and H is the Hubble parameter. We will only consider scattering with the Higgs bosons, as this 
will prove sufficient to demonstrate thermalisation. (Zero-momentum s particles could also annihilate with thermal 
background s with a similar rate.) The scattered s particles will achieve a thermal energy by subsequent scatterings 
and will achieve a thermal equilibrium number density via annihilation to SM particles, provided scattering occurs at 
a temperature greater than the freeze-out temperature of the s scalars. 

The scattering rate of zero- momentum s scalars by Higgs bosons is F « nav, where n is the number density of 
thermal background Higgs bosons, n = ^gnT^ with gn ~ 4, v ^ 1 is the velocity of the relativistic Higgs particles 
and a is the scattering cross-section, 



>2 



(30) 



Assuming ptotai — Prad, H ~ 3.4T^/Mp. Therefore thermalisation of the background will occur if 

^Ik > ,31) 



247r2 TUs^ Mp 



This is satisfied if 



m. 



< 1.0 X lO^^A^Mp . (32) 



This is easily satisfied, as we expect m^^ to be less than a few TeV. This means that any inflaton energy density 
remaining in the s oscillations will be scattered and thermalised as soon as the oscillation amplitude enters the 
oscillation regime, if not before. A key point is that since the energy density of the Universe is dominated by 
radiation at the onset of oscillations, there is sufficient energy to completely thermalise the oscillations. Of 
course, it is likely that thermalisation of the oscillations will occur even earlier; our calculation is simply a check that 
thermalisation of the s oscillations can safely occur. Thus there is no danger of a residual dark matter density in 
oscillations; s dark matter is entirely from thermal freeze-out of s particles. 



V. DISTINGUISHING S-INFLATION FROM HIGGS-INFLATION THROUGH OBSERVATIONS 

We next present the main result of our paper, which is the prediction of the spectral index n and tensor-to-scalar 
ratio r as a function of rrih for the three variants of the model. Compared to pure Higgs inflation, the gauge singlet 
variants of the model introduce two new parameters As and Xhs- The second can potentially be fixed by experiment. 
There will also be a minimum value of Xhs below which the s particles cannot account for the observed density of 
dark matter. We examine all values of As and Xhs sX ^ = rrit where the magnitude of all couplings (except ^0) remain 
perturbative (we generously take this to be < 100) up to the scale of inflation. We have not imposed any limits on 
the perturbativity of the potential in the non-inflationary direction. Along with stronger limits in the inflationary 
direction, this would impose stronger bounds on the permitted combinations of mh , As and A;is ■ Our aim here is only 
to compare the behaviour of n and r as a function of ruh ■ 

We flnd that the spectral index is well approximated (in most of the parameter space) by considering only the 
tree-level contribution to 77 



1 + 2^^1-^^. (33) 



Thus, the shape of n is determined by C^'/'^, where this is determined by the radiatively-corrected e through the 
normalisation to the COBE data (Eq. (|15p ) and through the integration to get TV, Eq. (|18p . If radiative corrections 
cause € to increase above its classical value, we would expect to be larger for a flxed N . This means that the 
magnitude of fj is increased (i.e. becomes less negative) and so n is increased. 

We will use the quantity as introduced in Section IHIl which is approximately the part oi due to radiative 
corrections, to explain how n deviates from its classical value. Large corresponds to large radiative corrections to 
e, which corresponds to larger n, as described above. Lh and Ls are given by Eq. (|2T|) and Eq. (|22|) . 
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(a) Spectral index n. The points curving downwards 
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FIG. 2: Spectral index n and tensor to scalar ratio r versus Higgs mass for inflation in the S-direction (pink circles), for inflation in 
the H-direction (blue circles) and for pure Higgs inflation (solid green line). Couplings have been varied by increments of 0.1 and 

required to remain less than 100. 



Spectral Index as a function of 



We first present our main result, which is the very different behaviour of the spectral index and tensor-to-scalar ratio 
as a function of Higgs mass in S'-inflation as compared with Higgs Inflation models. We then explain the behaviour 
of the spectral index as a function of the Higgs mass and model couplings in the different models. 

We show in Fig. (Pa|) the approximate range of n as a function of rrih for each model. In these plots As and Xhs are 
allowed to take any values that are multiples of 0.1. Introducing perturbativity and vacuum stability bounds would 
cause these areas to somewhat decrease in size. We see that there is a striking difference between the two models at 
larger values of m/j. At mh ^ 150 GeV, the possible values of n are very different^. There is more overlap at lower 
rrih — for 125 GeV ^ ruh ^ 135 GeV, it appears unlikely that n could discriminate between the models. Of course, 
a measurement of A^^ will leave only one free parameter (As), reducing the number of points available. 

We show in Fig. (j2bp the equivalent figure for r, with the same restrictions on As and Xhs as above. We see that r is 
in general low (r ^ 0.02) — although for large mh (s-direction) and small mh (i/-direction), it can take values which 
are only just within the current WMAP limit r < 0.22. Thus there is a small chance that r may be detectable by 
Planck, although these extreme points may be excluded when full stability and perturbativity constraints are applied 
to the models^". 

As we have shown, the error on the theoretical value of n, due to uncertainty in the reheating temperature, is 
An — ±0.001 or less. We expect the Planck experiment to measure n to a 2-a accuracy of ±0.0005. Therefore if 
Planck should find n significantly larger than 0.967 while LHC finds a Higgs with mass larger than 135 GeV, S'- 
inflation will be compatible with the observations but Higgs inflation will be ruled out. If the spectral index n is 
measured to be signiflcantly less than 0.965 then Higgs inflation will be compatible with observations while S'-inflation 
would be ruled out. 

Thus a combination of LHG data and Planck data can determine which weak scale inflation scenario is viable. The 
possible measurement of Xhs in direct detection dark matter experiments (and, if TOs < m/i/2, at the LHC), combined 
with a Planck measurement of n larger than (but close to) 0.967 and negligible r, would be strong evidence in support 
of S-inflation. 

We next explain in more detail the behaviour of the spectral index as a function of mh in the cases where the Higgs 
is the inflaton and s is the inflaton. 



The small number of s-direction points with n < n^i are likely to be ruled out by perturbativity and stability constraints on the 
couplings. 

These constraints should be applied in both directions {h and s), regardless of the direction of inflation, and should be applied at least 
to a scale just above the scale of inflation. 
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(b) Inflation in the Higgs direction, with As = 0.0 (solid 
red) and As = 0.15 (green dashed). 
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(c) Inflation in the S-direction with A/js = 0.0: As = 0.01 
(solid red) and As = 0.25 (green dashed). 'Pure' Higgs 
inflation is shown for comparison (solid black). 




110 120 130 140 150 160 170 180 
m/i ( GeV) 
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FIG. 3: Spectral index n versus Higgs mass for different versions of the model 



Spectral index with the Higgs as inflaton 

Firstly, we consider the effect of Xhs on Higgs inflation and so set Xs{mt) — (but allow for its running). The 
results are shown in Fig. ([5a| for Xhs = (red; pure Higgs inflation), A^^ = 0.1 (green), 0.3 (pink) and 0.5 (blue). 
The shape of the curves and range of n for all X^s are quite similar (for ruh > 122 GeV), but two main features can 
be seen: (i) as Xhs increases, the curves shift to the left, also shifting the range of mh and (ii) there is a turnover at 
low TO/i for Xhs = 0.3. These features are explained below. 

(i) A larger Xhs increases f3\^ (Eq. (jA-81) ). giving a larger Xh for a given nih- At large nih, Lh oc — 6A;i (Eq. ([HI)), 
and so n falls faster with larger Xhs- This explains the shift to the left as Xhs increases. The range of mh accessible 
to this model is shifted downwards as Xhs increases. A larger A^^ can prevent A^ becoming negative at low values of 
rUh- Large Xhs causes the couplings to grow faster and the perturbativity limit to be reached at a lower m/j. Thus, 
the range of mh is shifted. 

(ii) The turnover at low mh (seen for the case Xhs ~ 0.3 in Fig. (j3al) ) is due to the term in Lji oc Small Xh means 
that ^/j is small (from normalisation, Eq. (HJ). Large Xhs gives a large running of (remember that S,s(jnt) — 0) so 
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the ratio |^ is of 0(1) at the scale of inflation^^. Thus there is a large negative contribution (last term of Eq. (|2ip ) 
ahxiost balancing the positive term oc 

As shown in Fig. pb[) . As has a steepening effect on n. However, this is a relatively small effect for \hs{mt) — 0.3. 
Above about — 0.15, we reach the perturbativity limit of As- The steepening is due to As causing \hs to increase, 
therefore exaggerating the effects of Fig. pap further. The last term in Eq. ([?!]) . —2Xhs^, rnay also play a role. 
Increasing As increases ^s, which may give some contribution to the decrease in n at larger ruh- For Ahs('Tit) = 0.1 
(not shown) we find that increasing As has a negligible effect on n (As becomes non-perturbative at As(mt) ~ 0.3 — 0.4 
in this case). 

Wc conclude that introducing a real singlet scalar to the model of Higgs inflation can affect the spectral index 
prediction, increasing it at low m/i and decreasing it at higher m/j. This effect is controlled mainly by the magnitude 
of Xhs which is, in principle, measurable. The addition of A^s also changes the range of rrih, decreasing both upper 
and lower limits. Negative values of Xhs are allowed and give similar results. The spectral index is insensitive to As- 

Spectral index in S-inflation 

With Xhsi'mt) = we see that the spectral index does not vary noticeably with nih — see Fig. ([5c)) (where we show 
the pure Higgs inflation case for comparison). This is as we would expect, since with Xhs = 0, the model is completely 
decoupled from the Higgs sector. We also see that n does not vary with As (if As is increased much further than 
shown in the figure, it reaches its perturbativity limit). This is because the deviation of n from its classical value is 
determined by (Eq. which is proportional to A^^ = 0. 

With non-zero A/is, we show the results for Xs{mt) = 0.1 (solid) and Xs{mt) = 0.01 (dashed) in Fig. (pd)) . We see 

that increasing A^s has a dramatic effect on n. This is because cx -j^. Large \Xhs\ (at infiation scale) will therefore 
cause Ls to increase, giving larger n. Smaller As gives larger Ls, increasing n further, although the effect of Xhs is 
dominant. 

There is a significant sensitivity to the value of As. In S-infiation it is possible, in principle, to measure all the 
unknown couplings with the exception of As. Thus, although the qualitative behaviour of n as a function of m/i can 
be known, there is a fundamental obstacle to making a precise prediction of n for a given m/j in S'-infiation. However, 
with a reasonable estimate of the value of As, such as As ~ A;i, we can estimate the expected value of n. In addition, 
constraints on As from vacuum stability and perturbativity of the potential will also constrain the allowed range of n. 

VI. CONCLUSIONS 

The Standard Model and its minimal extension to include a real gauge singlet scalar dark matter particle can 
account for inflation if the scalar fields are non- minimally coupled to gravity with a large dimensionless coupling 
In this paper we have considered the precision predictions of these models for n and r as a function of Higgs mass. In 
order to ensure a consistent calculation and so avoid the contradictions found in the existing results for the spectral 
index in Higgs Inflation, we have recomputed the RG improved potential for Higgs Inflation and S'-inflation using 
the same method. The addition of the gauge singlet scalar allows inflation along both the Higgs direction (Higgs 
Inflation) and the singlet direction (S'-inflation). We have shown that these alternatives can be clearly distinguished 
by the deviation of the spectral index from its classical value as a function of mh if the Higgs mass is sufficiently large, 
ruh ^ 135 GeV. In S'-infiation, quantum corrections to the inflaton potential cause n to increase with nih, whereas in 
Higgs Infiation (both in the SM and in its singlet scalar extension), n decreases with increasing nih- The theoretical 
prediction of n is accurate to ±0.001, so observation of n by Planck, combined with an LHC determination of rrih, 
has the potential to distinguish which of the models is consistent with data. For smaller mh, both Higgs Infiation 
and S-inflation can produce values of n which are larger than the classical value, so in this case it is more difficult 
to distinguish the models, although since pure Higgs Infiation can, in principle, make arbitrarily accurate predictions 
for 71 as a function of nih, it would still be possible to exclude pure Higgs Inflation. We find that the additional 
singlet scalar also has a small but potentially observable effect on the predictions for n and r in the Higgs direction. 
Finally, all the models we have considered predict a small value of the tensor to scalar ratio r < 0.02 over most of the 
parameter space. 



However, with and ^h of similar magnitudes, we can no longer assume that inflation is in the _ff-direction. Therefore this is perhaps 
not a valid region of the parameter space — however we do not impose any constraint on 4^ in this chapter. 
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In most inflation models, the main theoretical uncertainty in the prediction of n is due to the reheating temperature. 
A striking feature of Higgs Inflation and its singlet variants is that the reheating temperature can be determined very 
precisely, since the relevant model couplings are either completely determined experimentally (in the case of pure 
Higgs Inflation) or determined up to the unknown values of and X^s (in S*- inflation) . Although there is a somewhat 
greater uncertainty in the case of S-inflation, with plausible assumptions based on the requirement of thermal relic s 
dark matter, a very small range of Tr can be determined. We flnd that n can be predicted up to an error An — ±0.001. 
Therefore if Planck can determine n to an accuracy of ±0.0005, then it should be possible to distinguish whether 
Higgs Inflation or 5*- inflation is consistent with observation for a range of mh- 

In the case of Higgs inflation, a precise test of the model is possible since n can be exactly calculated as a function 
of nih ■ In the case of S'-inflation, the s self-coupling probably cannot be directly measured but it can significantly 
modify the value of n. Therefore a precise prediction of n as a function of nih may not be possible. However, the 
magnitude of the positive deviation from the classical value of n can be estimated under plausible assumptions regard- 
ing the magnitude of As, including full consideration of the perturbativity constraints. Combined with observation 
of s dark matter in direct detection experiments, and possible confirmation of s particles by LHC in the case where 
TOs < mh/'2, this would provide strong support for S'-inflation. 

Our results are dependent upon a number of assumptions being correct. Firstly, we assume that the effective 
potential is not modifled by the physics that unitarizes Higgs scattering due to the non-minimal coupling to R. This 
is most likely to be true if the scale of perturbative unitarity violation, E A = Mp/^^, is in fact the scale at which 
strong coupling unitarizes the scattering process. In this case, no modification of the action is necessary. Secondly, in 
keeping with all previous analyses of the effective potential, we assume that the effect of the non-minimal coupling on 
the RG improved effective potential is fully taken into account by the suppression of the inflaton propagator in the 
Jordan frame at > A. The validity of both of these assumptions merits further investigation. 

Higgs Inflation and its variants generally assume that corrections due to Planck scale-suppressed operators, such as 
may be expected from the UV completion of the theory, do not signiflcantly modify the predictions of the model. Such 
corrections could modify both the inflaton potential and the non-minimal coupling. Whether these corrections can 
modify the model predictions for n and r will depend on the form of the UV completion. If the non-renormalizable 
potential terms in the Lagrangian produce physical interactions with strength determined by the Planck scale, then 
the corresponding Lagrangian operators will have additional factorial suppression factors which will ensure that their 
contribution to the spectral index is negligible. (Such corrections can be comparable to 1-loop corrections in the 
absence of factorial suppression [9].) Planck corrections to the non-minimal coupling will also be negligible if the 
full UV completion of the non- minimal coupling is an expansion of the form ^/(|iJp/Mp)i?, where / is expanded 
in |iJp/Mp. Since the effect of Planck-scale corrections on the effective potential is difficult to determine, the best 
strategy is to test the models experimentally. 

In conclusion, it is possible that inflation and thermal relic WIMP dark matter might be understood entirely in 
terms of the Standard Model and its simplest extension to include a singlet scalar dark matter particle. We have 
shown that the nature of inflation in such models may be determined using forthcoming LHC and Planck data for a 
range of Higgs mass and spectral index. 
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Appendix A: Radiative corrections and slow roll parameters 

In this appendix, we detail the equations firstly of the RG improved effective action (including the suppression of 
the scalar commutator) and secondly of the slow roll parameters. 

(i) The modifled commutator 

In our calculation of the RG-improved effective potential we have applied the prescription that the effect of the 
non-minimal coupling is to suppress the inflaton propagator in the Jordan frame when calculating Feynman rules. 
We briefly review the argument for this method. 
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The non-minimal coupling to gravity in the Jordan frame means that the scalar field propagator (proportional to 
the commutator [(j){x) , (j){y)]) is modified. The technique for calculating this modification was introduced by [33i] and 
first applied to the case of Higgs inflation by 14). 

The effective potential is calculated in the Jordan frame, using the RG equations derived in the Jordan frame, 
where the fields are canonically normalised. The scalar field commutator 

(A-1) 

is computed by first transforming the action to the Einstein frame, where the gravitational term is minimal and the 
canonical momentum tt can be calculated. This is given by 

^ V V f S V (5^^5.0) . (A-2) 



where x is the rescaled field with canonical normalisation in the Einstein frame. Using the commutator Eq. (jA-l|) 
and rearranging, the result is 

= ihc4,S'-^\x-y). (A-3) 

Therefore all propagators for scalars (jn ^re in principle suppressed by factors c^^ = (^f^) ' However, the only 
scalar propagator which is suppressed is the non-minimally coupled inflaton field which has a large expectation value, 
since for all other scalar fields (^^^ = jp-- (This can be seen from Eq. ([5]).) In this case the commutator should be 
understood as applying to perturbations about the background infiaton field. The suppression factor for the inflaton 



is given by [33|, \34 



1 I 



c<* = ^-T^- (A-4) 

l + (6e, + l)^ 

It is important to emphasise that in the case of the Higgs as inflaton, only the physical Higgs scalar h has a suppression 
factor, not the Goldstone bosons in H. In our calculations, although suppression factors for both h and s are included 
in the RG equations, the suppression factor for the field which is not the inflaton is set to 1. 



(ii) RG equations 

We use the two-loop RG equations for all the SM couplings, inserting the suppression factor for each non- 
minimally coupled inflaton running in a loop, as described in the previous section. For any coupling A we have 

'^^f. (A-5, 

for inflation in the s direction (since 7s = 1), and 

dX Px 



dt l + jH 



for inflation in the h direction, where t = In-^ and 



1 ^9g2 3.g'2 ^ ,\ 1 /271 4 9 , ,^ 431 ,4 



(A-6) 



= -m^[— + — -^y')-Jl6^[-32' "16^^ "W^"^ 

-I {Ig' + + 8-93) + T^'^yt " 64A^.) . (A-7) 
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The SM one- and two- loop equations can be found in and [2g| . The RG equations for the scalar couplings can 
be obtained""^^ using the technique detailed in 35[, as in Q- The one-loop /3-functions for the scalar couplings are 



-6yt + I {2g' + {g^ + 5'^)^) + {-9g' - Zg'^ + Uy^) A, + {iScl + 6) Xl + ^c^AL (A- 



167r2/3W = 4chC,AL + 6 (c^ + l) A,,A,,, - - {3g^ + g'^) + Qy^Xhs + GclXsXhs 



(A-9) 



and 



The RG equations for the non- minimal coupling can be derived following [s^, as in 0] (see also 



2 \2 



(A-10) 



16^'^ = (3 + Ch) Xhs Uh + + + 



(A-11) 



and 



167r2% = ( (6 + 6c,0 A„ + 6y^ - ^{3g^ + 5") ) ( ^ + ^ ) + ( 6 + ^ ) c.A„, 



dt 



(A-12) 



The equations for the gauge and Yukawa couplings are 

.3 /3 



/3. 



39-Ch 3 g 
-9 



12 



(A-13) 



12 ^ ' 167r2 V 18 ' 2 ' 3 6 / 



(A-14) 



and 



yt 



9 
i 

1187 



5 



17 
12 



23 ^ 2 \ 2 V 2/t 



23 4 3 2/2 

- iff 



^9" + Q9^9l + Ys'al - 108ff^ + C^g^ + ^5" + 3633^ ) c^y^ 



+6 {-2clyt - 2clyfXh + clXl) 



(A-15) 



(iii) Slow roll parameters 

The Jordan frame RG improved effective potential of Eq. ^ is transformed to the Einstein frame in order to 
compute the slow-roll parameters. Its dependence on fi = (j) during inflation is then given by the variation of the 
couplings from RG equations as well as the explicit dependence on the inflaton (j). We can then relate the slow-roll 
parameters of the quantum corrected effective potential directly to the RG equations. 

Starting with the definition 



( 1 dU 



U dx4, 



Ml /d^ 
2 VrfX0 



^ 'dU 1 ^ ^ 



(A-16) 



These were also obtained by although there is a minor difference between our work and Q relating to a factor of Cf^ in the RG 
equation. The term in proportional to is suppressed by a factor (1 -|- cjj) in Eq. (A2) of but (1 -|- c^) in our work Q. We 
believe that the latter is correct because only one physical Higgs h runs in the corresponding loop. In practice the effect of the difference 
between these is negligible. 
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u d(j) ^ (jiff^ (j) yx^ dt n^{i + -f^) n^M^ dt J 

to give 

_ ^ Ml fd^y ( ^ ]_dX^ _ 47^ _ 2£Gld^' 

" 2(j)^\dx^) \^ dt n^i + j^) n^M^ dt ^ 

It should be noted that 7^ = and G, = 1. We have also used ^ = - ^i+ll"t)^ ■ 
To find 77, we start with the definition 

^ _ M2 d^u _ M2 fd(t>y d^U d'^(j) dU 



(A-17) 



2 



(A-18) 



U dxl U \dx) d(j)^ U dx^ dcf) 



(A-19) 



As this expression will be complicated, we assume = is a good approximation during inflation (when we evaluate 
n and r), since ~ from Eq. (IA-4[) . Then we have 

I d^U _ lldU 1 (duV Se^G^ 8G2 dS,^ 



U d(j? (f)U d<f) \d<f) J n^M^ (1 + "fh)^^M^ dt 

2(/)2g4 /d^^V 2G2 d^^^ 8eG2(27,, + 7,2) 



n^M^ \dt J n^M^ dt^ mM2(l+7;,)2 
1 / 1 fdXhV , 1 d^Xh 4 ^7,,' 



'2 I Xl\dt ) Xh df^ 9?{l+-ihf dt 



(A-20) 



where 



d^Xh ^_drff^dXh_ 1 

df^ ~ 1 + 7,, dt dt 167r2(l+7,0 



„ , dXh „ , r, dyt 9 q do 



2^ 2-^^ dt 2-^ ^ dt ^ ^ ^ dt 



(A-21) 



and 



drfh ^ 1 / 9g dg 3g' dg' 
dt ~ 167r2 2 dt 2 dt 
9 

dt 



6y, 



dt 



1 



2 ,dg' 5 I 9 dg 17 ,dg' 



2 V 2"^ dt 



6'' dt 

dyt ' 
dt 



16g: 



(167r2)2 
dg3 



271 ^dg 9 ,2dg 



dt 



dt 



2/? 



99' 



dt 



(A-22) 



dt2 + 1/6 V dt 



^0 + 1/6 
167r2 



g rfAfe ^ i2yf— - 9g— - 3g' — 
dt dt dt ' dt 



(A-23) 



In these expressions, 



d20 _ 1 /d0\" dr? 
dP ~ (v^. 



1 + 3A/2 4 



dl^V d2(l]2^ 



d^j 



(A-24) 
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where 

dn 1 ^(j)G^ / 27 1 d(, 



2n ivn \ 1 + 7 c dt 



and 



l + 7/i 




(A-25) 



(A-26) 



Appendix B : Reheating in ^-inflation 

Reheating in S'-inflation is similar to reheating in Higgs Inflation. In both models reheating occurs primarily through 
preheating via stochastic resonance of the inflaton to bosons, which subsequently annihilate to relativistic Standard 
Model fields. In the case of Higgs Inflation this occurs via excitation of gauge boson modes by the oscillating Higgs 
inflaton, Xh, in which case the coupling is the (known) weak gauge coupling g. In the case of S'-infiation the oscillating 
Xs field excites the Higgs bosons via the (unknown) coupling A/is. In the analysis of 0, the gauge boson modes are 
approximated by scalar modes. This analysis may therefore be directly applied to the case of Higgs boson modes 
excited by the oscillating Xs field. It is also possible for the infiaton to excite its own modes via the self coupling 
(A;i and As for Higgs Inflation and S'-inflation respectively, where A^ is known and A^ is a free parameter). In the 
following we will adapt the method discussed in the appendix of 9] to the case of S-inflation, discussing only the key 
results and the differences. 

In S'-inflation, Higgs bosons are created when Xs passes through zero at times tj. Between tj^i and tj the occupation 
number of mode k remains constant at nj,. This changes to n'j^^ after Xs passes through zero at tj, where Q 



< - ^tr ^ 'V 1 + <^<fBf.<^osiou (B-1) 

Here Rk and Dk are 'reflection' and 'transmission' coefficients relating coefficients and f3j, of the Higgs mode 
function = a^^'^4'k before and after Xs passes through zero at tj, and Q-j-^f = —291. ~2(i'^^ + f) + ^'"S'^i ~ ^-^g/^L 
where O^, is the phase of the mode function. There are two regimes, depending on the size of n^. 

(i) Uk <^ 1: During the time At = ^ between zero crossings, n;^ ~ jsTp ^'^^^^^^ 4'i are produced. (Here (f)i {i — 1,...4) 
arc the four real scalars in the Higgs doublet.) The majority of these particles decay (to relativistic Standard Model 
particles) before the next zero crossing and there is no significant transfer of energy to radiation. 

(ii) rifc 3> 1: it is now the second and third terms of Eq. ()B-1I) that are important and we can write (for each ipi) 

where fik is an average over the angle OlotJ which is assumed to be completely randomly distributed. This can be 
calculated, giving 



1 

2^ 

1 



In 



\Rk\ + i 
\Rk\ 



^ +i-\n{\Dk\^) 



^ln(l^fcP) . (B-3) 

In order to have a resonance, the reflection coefficient, must be reasonably large. This means that \Dk\ is small 
(from the requirement + |£'feP = 1) and this gives the final line of Eq. ()B-3p . We are interested in the rate of 

change of the total number of particles produced, so we must integrate over k. Thus, as the change in nk during the 
time between zero crossings of the infiation. At — ^, is Auk — 2'K^k'nk, the total rate of production of Higgs scalars 
is 

dn oj f°° , d^k 



18 



where n is the total number of (fii scalars. This integral is not calculated exactly, instead it is assumed that /Xfc is well 
approximated by fi{k = 0) = B Q, which is found to have the value B = 0.045. Summing all four real Higgs bosons 
in H gives for the rate of change of total number density of Higgs bosons tit 

— - 4 X 2ujBn = 2ujBnT. (B-5) 
dt 

In the components of the and Z bosons produced via preheating by the Higgs boson are treated as scalars. 
This is not strictly true for that case but it is correct for our calculation of the production of scalar Higgs bosons via 
preheating. 

The Higgs bosons which arc produced at the zero crossings can decay in the adiabatic regime between zero crossings 
where their mass is large. The Higgs has a large mass because the background Xs held is large. Therefore it can decay 
via the Yukawa coupling to fermion pairs, primarily the top quark. The decay width of the Higgs bosons (with mass 
averaged over an Xs oscillation) is given by [37l | 

327rm„/ ^ ^ 17.65 ^ 



where {mi) = ( ''^''^^ ). The Higgs bosons are non-relativistic when produced but their decay and annihilation products 



2 

are relativistic. 



Denoting the amplitude of the oscillation of x by X, we can investigate whether reheating can complete before 



quartic oscillations begin, defined to be Xqr = Y I ■ O^i^e the decay of the produced Higgs bosons becomes 
subdominant to their production, an exponential regime can begin, corresponding to Bose-enhanced generation of 
Higgs bosons. Up to this point, no significant energy transfer from Xs to the Higgs bosons occurs since Higgs decay 
prevents a large Higgs occupation number. Reheating completes fairly rapidly after this, as once the Higgs bosons 
are produced through stochastic resonance, they can annihilate and produce a thermal background, draining most 
of the energy from the Xs oscillations and causing its amplitude X to rapidly decrease to Xc r at which point the 
quartic oscillations begin. The upper limit on the time of reheating is the time at which the decay and the exponential 
production are equal. Using Eq. (|B-6|) and Eq. (jB-5|) we find that the two processes are equal at 

^ ^^Xc« 4.0^ XcH ■ (B-7) 

-^hs ■^hs 



This is the absolute maximum value of X at which reheating occurs via production of H bosons. After Xs — Xsto, 
Xs will rapidly decrease to Xcr., therefore prad > Pinf at Xcb- The calculation assumes instant annihilation of the 
_ff-bosons to relativistic particles and a very efficient stochastic resonance, therefore Xsto is an upper limit of the value 
of X at which reheating occurs. For this reason all estimates of Tr from stochastic resonance are upper limits. It is 
interesting to note that increasing the Higgs-s coupling, A^^, causes reheating to occur later. This is because a large 
\hs means that m| is large, so the Higgs bosons can decay for a longer period, preventing the exponential regime 
from beginning and therefore delaying reheating. 

The time of reheating is dependent on the ratio of couplings . Requiring this process of reheating to occur before 
X = XcR gives 

A, > 0.25A,,, . (B-8) 



If this condition is not fulfilled, reheating is not ruled out. It could either occur by directly producing excitations 
of the inflaton, as discussed below, or possibly during the quartic potential regime. Since our goal is to show that 
reheating in S-inflation can occur for natural values of the couplings, we will conservatively assume that Eq. (IB-Sp 
is satisfied. If it is not satisfied and if reheating can occur in the quartic regime, a somewhat lower temperature of 
radiation domination Tr would be expected. 

An alternative method of draining energy from the background field is the direct production of excitations of the 
inflaton Xs, which is possible as the inflaton potential changes at small Xs to V xt- The mode equation can be 
solved perturbatively, assuming the number of particles produced is small. The s particles produced are relativistic, 
with energy density ^ 

3 u}^ 

Pexcitation — TT 7^ ■ (-^~^) 
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From this, if the production of excitations is the only process of reheating then radiation domination occurs when the 
Xs oscillation amplitude is X = X^x , where 

Xex ^ 7y/KXcB.- (B-10) 

If this mechanism of reheating is to occur before the inflaton potential becomes quartic, we require 

As > 0.019. (B-11) 

The relativistic Xs particles are expected to subsequently annihilate to Higgs bosons, so producing the thermal 
background of SM particles. This is most efficient for large Xhs- Both of these limits are shown in Fig. ([T]). 
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